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Physics of current fluctuations has proven, during the last 15 years, to be a very profound topic of 
electron transport in mesoscopic conductors (for a review, see ref. [T]). Usually, current fluctuations are 
characterized by their spectral density S^w) measured at frequency to: 

S 2 {u>) = {i{u)i{-u>)), (1) 

where i(u>) is the Fourier component of the classical fluctuating current at frequency u> and the brackets 
(.} denote time averaging. In the limit where the current can be considered as carried by individual, 
uncorrelated electrons of charge e crossing the sample (as in a tunnel junction), S^w) is given by the 
Poisson value ^(w) — el and is independent of the measurement frequency u>. At sufficiently high 
frequencies, however, this relation breaks down and should reveal information about energy scales of the 
system. In particular, in the quantum regime fao > eV (V is the voltage across the conductor), it turns 
out that the noise cannot be seen as a charge counting statistics problem anymore even for a conductor 
without intrinsic energy scale. In this regime, the noise spectral density reduces to its equilibrium value 
determined, at zero temperature, by the zero-point fluctuations (ZPF): 

Si eq) (u;) = Gfiw, (2) 

with G the conductance of the system. Experimental investigations of the shot noise at finite frequency 
have clearly shown a constant (voltage independent) noise spectral density for fao > eV in several systems 
El Sj ■ Although these experiments were not able to give an absolute value of the equilibrium noise 
(because of intrinsic noise of linear amplifiers used for the measurement), one has good reasons be believe 
that ZPF can be observed with this kind of amplifiers. Indeed, it as been proven in other detection 
schemes, theoretically 0E] and experimentally [3 El El HO], that ZPF can be detected from deexcitation 
of an active detector whereas they cannot be detected by a passive detector which is itself effectively in 
the ground state. 

In view of recent interest in the theory of the full counting statistics (FCS) of charge transfer, attention 
has shifted from the conventional noise (the variance of the current fluctuations) to the study of the higher 
cumulants of current fluctuations. Whereas the discrimination between active and passive detector seems 
to be clear for noise spectral density measurement, the situation is more complex for the measurement 
of high order cumulants at finite frequency. Indeed, the issues of detection scheme are closely related to 
the problem of ordering quantum current operators and, if the problem can be solved in a general way 
for two operators [HI 111) , measurements of higher cumulants are pointing out the problem of appropriate 
symmetrization of the product of n current operators: 

S„(w) = (i(uoi) i(uj 2 ) ...i{oj n )) 8(u>i + uj 2 + ... + lu„) (3) 

It is the goal of this paper to clearly present the problem of the third cumulant measurement on a well 
defined experimental setup using a linear amplifier as a detector. Until now, measurements of the third 
cumulant S3 of voltage fluctuations have been performed at low frequency, i.e. in the classical regime 
fao < eV, ksT where voltage fluctuations arise from charge transfer process [THEI1H1]- ^ e re P ort nere 
the first measurement of S3 at high frequency, in the quantum regime fan > eV, ksT. It raises central 
questions of the statistics of quantum noise, in particular: 
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1. The electromagnetic environment of the sample has been proven to strongly influence the measure- 

ment, through the possible modulation of the noise of the sample [12] . What happens to this 
mechanism in the quantum regime? 

2. For eV < huj, the noise is due to ZPF and keeps its equilibrium value: 52 = Ghw with G the 

conductance of the sample. Therefore, S2 is independent of the bias voltage and no photon is 
emitted by the conductor. Is it possible, as suggested by some theories [13 [T()J [T7] , that 53 7^ in 
this regime? 

In the spirit of these questions, we give theoretical and experimental answers to the environmental 
effects showing that they involve dynamics of the quantum noise. We study the case of a tunnel junction, 
the simplest coherent conductor. Using these results, we investigate the question of the third cumulant 
of quantum noise. 



1 Environmental Effects and Dynamics of Quantum Noise 

We show in this section that the noise dynamics is a central concept in the understanding of environmen- 
tal effects on quantum transport. First, we present a simple approach (in the zero frequency limit) to 
calculate the effects of the environment on noise measurements in terms of the modification of probability 
distribution P(i) of current fluctuations. We do not provide a rigorous calculation, but simple consider- 
ations that bear the essential ingredients of the phenomenon. This allows us to introduce the concept of 
noise dynamics and determine the correct current correlator which describes it at any frequency. Then, 
we report the first measurement of the dynamics of quantum noise in a tunnel junction. We observe that 
the noise of the tunnel junction responds in phase with the ac excitation, but its response is not adiabatic, 
as obtained in the limit of slow excitation. Our data are in quantitative agreement with a calculation we 
have performed. 



I(t) = <l> + i(t) 
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Figure 1: Schematics of the experimental setup. Current fluctuations i(t) = I(t) — {I) are measured by 
an ampmeter with a bandwidth A/. 



1.1 Effects of the environment on the probability distribution P(i) 

In the zero-frequency limit, high order moments are simply given by the probability distribution of the 
current P(i) calculated from the current fluctuations measured in a certain bandwidth A/ (see figO}: 

M„ = Jt n P(t)di (4) 

The cumulant of order n, S n is then given by a linear combination of Mk A/ fc_1 , with k < n [18j . In 
practice, it is very hard to perfectly voltage-bias a sample at any frequency and one has to deal with the 
non-zero impedance of the environment Z (see Fig. [l}. If V fluctuates, the probability P(i) is modified. 
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Let us call P(i] V) the probability distribution of the current fluctuations around the dc current / when 
the sample is perfectly biased at voltage V, and P(i) the probability distribution in the presence of 
an environment. R is the resistance of the sample, taken to be independent of V. If the sample is 
biased by a voltage Vq through an impedance Z, the dc voltage across the sample is V — Rn/ZVo with 
i?H = RZ / (R + Z). The current fluctuations in the sample flowing through the external impedance induce 
voltage fluctuations across the sample, given by: 

/> + oo 
Z(uo)i(uo)e lu)t duo (5) 
-oo 



Consequently, the probability distribution of the fluctuations is modified. This can be taken into account 
if the fluctuations are slow enough that the distribution P(i) follows the voltage fluctuations. Under this 
assumption one has: 

dP 

P{i) = P(i; V + SV) ~ P(i, V) + SV — + ... (6) 

ov 

supposing that the fluctuations are small (SV <C V). One deduces the moments of the distribution (to 
first order in SV = —Zi) for a frequency independent Z: 

M n ^M n -Z d -^ + ... (7) 



This equation, derived in Ref. |19j . shows that environmental correction to the moment of order n is 
related to the next moment of the sample perfectly voltage biased. For n = 1 we recover the link between 
noise and Dynamical Coulomb Blockade through the noise susceptibility (see below) that appears as 
dM^jdV in the simple picture depicted here [2Dj. Let us now apply the previous relation to the third 
cumulant (S n = M n A/™^ 1 for n = 2,3): 

S 3 ^S 3 -3Z^cS 3 -3ZS 2 ^- r (8) 



It is a simplified version of the relation derived in refs. |21| . The way to understand this formula is the 
following: the first term on the right is the intrinsic cumulant; the second term comes from the sample 
current fluctuations i(i) inducing voltage fluctuations across itself. These modulate the sample noise £2 
by a quantity —Zi(t)dS2/dV . This modulation is in phase with the fluctuating current i(t), and gives 
rise to a contribution to the third order correlator (i 3 (t)). This environmental contribution involves the 
impedance of the environment and the dynamical response of the noise which, in the adiabatic limit 
considered here, is given by dS2/dV. However, at high enough frequencies, and in particular in the 
quantum regime Hoj > eV, this relation should be modified to include photo-assisted processes. The 
notion of dynamical response of the noise is extended in the following section to the quantum regime in 
order to subtract properly the environmental terms in the measurement of the third cumulant. 



1.2 Dynamics of Quantum Noise in a Tunnel Junction under ac Excitation 

In the same way as the complex ac conductance G(ujq) of a system measures the dynamical response of 
the average current to a small voltage excitation at frequency loq, we investigate the dynamical response of 
current fluctuations x^ (w), that we name noise susceptibility. It measures the in-phase and out-of-phase 
oscillations at frequency ujq of the current noise spectral density ^(w) measured at frequency uj. In order 
to introduce the correlator that describes the noise dynamics, we start with those which describe noise 
and photo-assisted noise. Beside the theoretical expressions, we present the corresponding measurements 
on a tunnel junction [22] . It allows to calibrate the experimental setup and give quantitative comparisons 
between experiment and theory. 

Noise and photo- assisted noise 

The spectral density of the current fluctuations at frequency w of a tunnel junction (i.e. with no internal 
dynamics) biased at a dc voltage V is [1] : 

ft(^)= s ° ( " +) : 5S( "- > . (9) 
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Figure 2: Top: Measured noise temperature T/v = S2{to) / (AksG) of the sample plus the amplifier with 
no ac excitation. Bottom: measured differential noise spectral density dS2(to)/dV for various levels of 
excitation z — e8V/(Htoo). z ^ corresponds to photo- assisted noise. Solid lines are fits with Eq. (fTU|) . 



where to± = to±eV/h. S^ito) is the Johnson-Nyquist equilibrium noise, Sv>( w ) = 2G7kJ coth (hto / (2ksT)) 
and G is the conductance. At low temperature, the 5*2 vs. V curve (obtained at point C on Fig. [3]) has 
kinks at eV = ±fiw, as clearly demonstrated in our measurement, see Fig. [5] top. The temperature of 
the electrons is obtained by fitting the data with Eq. [9l We obtain T = 35 mK, so that Huj/ksT ~ 8.5. 
Note that a huge, voltage independent, contribution T/v ~ 67 K is added to the voltage dependent noise 
coming from the sample which masks the contribution from the ZPF. When an ac bias voltage SV coswoi 
is superimposed on the dc one, the electrons wavefunctions acquire an extra factor J2 n J n (z) exp(intoot) 
where J n is the Bessel function of the first kind and z = eSV/ (htoo) . The noise at frequency lo is modified 
by the ac bias, to give: 
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Jl(z)S 2 {V - nhuo/e, lo) 



(10) 



This effect, called photo-assisted noise, has been measured for to — 2 . We show below the first 
measurement of photo-assisted noise at finite frequency to. The multiple steps separated by eV = Hloq 
are well pronounced and a fit with Eq. [10] provides the value of the rf coupling between the excitation 
line and the sample SV (see Fig. [2] bottom). 

Noise susceptibility 

Photo-assisted noise corresponds to the noise S^(u) in the presence of an excitation at frequency too, 
obtained by time averaging the square of the current filtered around to, as in [5J for to = and in 
[3] for to ~ wo- This is similar to the photo- voltaic effect for the dc current. The equivalent of the 
dynamical response of current at arbitrary frequencies loq is the dynamical response of noise at frequency 
ujq. It involves the beating of two Fourier components of the current separated by itoo expressed by 
the correlator (i(v)i(u>Q — to)) . Using the techniques described in [T], we have calculated the correlator 
that corresponds to our experimental setup, using the "usual rules" of symmetrization for a two current 
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Figure 3: Experimental setup for the measurement of the noise dynamics X(loq, u>) and the third cumulant 
63(0;, wo — to) for u> ~ ujq. The symbol © represents a combiner, which output is the sum of its two inputs. 
The symbol (g) represents a multiplier, which output is the product of its two inputs. The diode symbol 
represents a square law detector, which output is proportional to the low frequency part of the square of 
its input. 



correlator and a classical detector. We find the dynamical response of noise for a tunnel junction [2D] : 

X(u 0) uj) = ^ J~] J n (z)J n+ i(z) (S$((J+ + nuj ) - S%(uj- - nui )) (11) 

n 

Note the similarity with the expression giving the photo-assisted noise, Eq. (TTT)|) . Note however that 
the sum in Eq. (|11[) expresses the interference of the processes where n photons are absorbed and nil 
emitted (or vice- versa), each absorption / emission process being weighted by an amplitude J n (z) J n ±i(z). 

Experimental setup 

The sample is an Al/Al oxide/Al tunnel junction identical to that used for noise thermometry [55]. We 
apply a 0.1 T perpendicular magnetic field to turn the Al normal. The junction is mounted on a rf sample 
holder placed on the mixing chamber of a dilution refrigerator. The resistance of the sample Rq = 44.2 ft 
is close to 50 ft to provide a good matching to the coaxial cable and avoid reflection of the ac excitation. 
The sample is dc voltage biased, ac biased at ujq/2it — 6.2 GHz, and ac coupled to a microwave 0.01-8 
GHz cryogenic amplifier. To preselect the high-frequency component i(uS), we use a 5.7 — 6.7 GHz band- 
pass filter (Fig. [3l lower arm). Its beating frequency uj is shifted to low frequency Su> by using a square 
law detector and the reference signal Vq cos(uj(jt + (p) in order to mix it with the low-frequency component 
i(Su). The power detector has an output bandwidth of Stu/2ir ~ 200 MHz, which limits the frequencies 
to contributing to the signal: \u>\ G [luq — Suj, ujq + Slu]. The low frequency part of the current, at frequency 
ui — u)q, is selected by a 200 MHz low pass filter (Fig. 03 upper arm). 

Experimental results 

We could not determine the absolute phase between the detected signal and the excitation voltage at the 
sample level. However we have varied the phase if to measure the two quadratures of the signal. We 
have always found that all the signal can be put on one quadrature only (independent of dc and ac bias, 
see inset of Fig. 0] (b)), in agreement with the prediction. In the case of as small voltage excitation, we 
define the noise susceptibility which is for noise the equivalent of the ac conductance for current: 

Xluq ( w ) expresses the effect, to first order in SV, of a small excitation at frequency u>o to the noise measured 
at frequency u>. We show in Fig. H] the data for X(ujq,uj)/SV at small injected powers as well as the 
theoretical curve for Xu> a (u = u>q): 

Xfc ,M = X *(0) = JL (£>+) - S>_)) (13) 
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All the data fall on the same curve, as predicted, and are very well fitted by the theory. The cross-over 
occurs now for eV ~ faj. Xu{u) is clearly different from the adiabatic response of noise dS2(co)/dV 
(solid line in Fig. HJ. However, in the limit SV — > and uio — > (with z <C 1), Eq. (|13[) reduces to 
Xw(0) ~ (l/2)(dS2/dV). The factor 1/2 comes from the fact that the sum of frequencies, ±(w + wo) 
(here ~ 12 GHz), is not detected in our setup. This is the central result of our work: the quantum noise 
responds in phase but non-adiabatically. 
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Figure 4: Normalized noise susceptibility Xu>{^) v $- normalized dc bias. Symbols: data for various levels 
of excitation (z = 0.85, 0.6 and 0.42). Dotted and dashed lines: fits of x^{ u ) (Eq- (fT5|) ). Solid line: 
(l/2)dS2/dV (experimental), as a comparison. Inset: Nyquist representation of X(ujq,uj) for z = 1.7 (in 
arbitrary units). The in-phase and out-of-phase responses are measured by shifting the phase ip of the 
reference signal by 90°. 



2 Third Cumulant of Quantum Noise Fluctuations 

2.1 Operator odering 

A theoretical framework to analyze FCS was developed in Ref. [23] to evaluate any cumulant of the current 
operator in the zero- frequency limit. In order to analyze frequency dispersion of current fluctuations it 
is necessary to go beyond the usual FCS theory [T5l IT6| [17] . An essential problem in these approaches 
is to know what ordering of current operators i corresponds to a given detection scheme. This problem 
is simpler for S2: the correlator S+(u>) — (i{uj)i{~uj)) with lj > represents what is measured by 
a detector that absorbs the photons emitted by the sample, like a photo- multiplier. The correlator 
S-(lu) — u>)i(u>)) = S+(—lu) represents what the sample absorbs, and can be detected by a detector 
in an excited state that decays by emitting photons into the sample. Finally a classical detector cannot 
separate emission from absorptions, and measures the symmetrized quantity: 

s s 2 ym -(uj) = (IMiiz^ll ± ('(-^M) ( 14 ) 

However, according to Kubo formula S+(w) — S-(ui) — Ghio , all these contributions have similar voltage 
and temperature dependence, at least for a linear conductor. In contrast, different ordering of three 
current operators give rise to very different results. The prediction for the Keldysh ordering, which is 
supposed to correspond to a classical detection, is: £3(0;, u>') = e 2 /, independent of frequency even in the 
quantum regime. As far as we know there is no clear interpretation of this ordering in terms of absorption 
and emission of photons. We give below two detection schemes for the measurement of S^w, 0) that may 
lead to different results. 
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2.2 Measurement of S v s(u,0) 



Experimental setup 

We have investigated the third cumulant S v 3(lu,0) of the voltage fluctuations of a tunnel junction in 
the quantum regime hw > eV. For technical reasons (the input impedance of the rf amplifier is fixed at 
Z — 50 fi), we measured voltage fluctuations v(t) instead of current fluctuations i(t). Thus, the impedance 
i?H = RZ/(R + Z) will act as the environment and will affect the measurement of the third cumulant 
83(10, 0) of the current fluctuations. We use the same experimental setup and sample as for the noise 
dynamics measurement, the only change is that the ac excitation is switched off: SV = Vo = (see Fig. 
[51 Thus only the noise of the amplifier can modulate the noise of the sample. A 5.7 — 6.7 GHz band-pass 
filter followed by a square law detector allows to mix high-frequency components v(uj) v(—lu — 5lo) which 
are multiplied by low-frequency components selected by a 200 MHz low pass filter, we end up with a dc 
signal proportional to SV 3 °c ( W ( LJ ) v{—uj — Suj) v(Sui)). The fact that the same setup is used to detected 
S3 and x is quite remarkable: it clearly indicates that the environmental effects in S3 are indeed described 
by x and not by dS2/dV. 




Figure 5: (a) Measurement of S v 3(uj,0) vs. bias voltage V (circles). The solid line corresponds to the 
best fit with Eq. (fT5|) . The dash dotted line corresponds to the perfect bias voltage contribution and the 
dotted lines to the effect of the environment, (b) Measurement of S3 (a;, 0) vs. bias voltage V (squares). 
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Experimental results 

S v s at T = 35 mK is shown in Fig. [5] (a), these data were averaged for 4 days. These results are clearly 
different from the voltage bias result because of the environmental contributions. As described before 
(see section 1), the noise of the sample is modulated by its own noise and by the noise of the amplifier 
S 2 ,n, to give rise to an extra contribution to S v s. By generalizing the expression (J5J, we find, assuming 
real, frequency independent impedances to simplify the expression (but we used the full expression for 
the fits of the data): 

S>(^,0) = -iZffikM) + Rj (52,jv(0) + 5 2 (0))xo(w) + 

(15) 

+ i?f (VM+5 2 H)x,(0) + iZf {S 2 ,n(") + xU") 

To properly extract the environmental effects, we fit the data obtained at different temperatures (35 mK, 
250 mK, 500 mK, IK, 4.2 K). The parameters R\\(0), R\\(u), S 2lN (0) and S 2 ,jv(0) that characterize the 
environment are independent of temperature, whereas S 2 (V) and x(Y) have temperature dependent 
shapes. This allows for a relatively reliable determination of the environmental contribution. We have 
performed independent measurements of these parameters and obtained a reasonable agreement with 
the values deduced from the fit. However more experiments are needed with another, more controlled 
environment, to confirm our result. The intrinsic S3 in the quantum regime, obtained after subtraction 
of the environmental contributions, is shown in Fig. [5] (b). It seems to confirm the theoretical prediction 
by [TEIEI (solide line), i.e. Sa((J, 0) = e 2 I even for fku > eV. 
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Figure 6: (a) Experimental detection scheme. The symbol (^) represents a multiplier, which output 
is the product of its two inputs. The diode symbol represents a square law detector, which output is 
proportional to the low frequency part of the square of its input. S^^uj, 5u> — » 0) is given by the product 
of the square of high frequency fluctuations with low frequency fluctuations, (b) Equivalent detection 
scheme using a photodetector to measure square of high frequency fluctuations. 



3 Conclusion 

We have shown the first measurement of the noise susceptibility, in a tunnel junction in the quantum 
regime Two ~ Jtuq 3> fc^T (uj/2ir ~ 6 GHz and T ~ 35 mK) 0]. We have observed that the noise responds 
in phase with the excitation, but not adiabatically. Our results are in very good, quantitative agreement 
with our prediction based on a new current-current correlator Xw ( w ) ^ (i(w)i(wo — w)). Using the fact 
that the environmental contributions to S3 are driven by Xj we have been able to extract the intrinsic 
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contribution from a measurement of (v 3 ) on a tunnel junction in the quantum regime. Our experimental 
setup is based on a "classical" detection scheme using linear amplifiers (see Fig. [6] (a)) and the results 
are in agreement with theoretical predictions: S^w, 0) remains proportional to the average current and 
is frequency independent [TH [T7j . This result raises the intriguing question of the possibility to measure 
a non-zero third cumulant in the quantum regime Hui > eV whereas the noise 62(0;) is the same as at 
equilibrium, and given by the zero-point fluctuations. 

One can think of another way to measure 63(0;, 0) with a photodetector (sensitive to photons emitted 
by the sample), as depicted in fig. [6] (b). In this case S3 is the result of correlations between the low 
frequency current fluctuations and the low frequency fluctuations of the flux of photons of frequency 
u> emitted by the sample. Since no photon of frequency u> is emitted for eV < fkv, the output of the 
photo-detector is zero and 83(01, 0) = 0. The expectation of such a measurement is sketched by a dashed 
line in fig. 03(b). Note that such a detection scheme has already been applied on laser diodes [2"4"1 I25j. 
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